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Abstract 

We study a blend of two kinds of homopolymers with tendency for seg- 

ragation. Cross-links between the chains of different kinds do not allow 

macrophase separation. Instead microphase structure appears. Starting from 

a microscopic model we derive the effective Hamiltonian and calculate the 

form of the inverse scattering function and the domain size. The latter is 

found to be of the order of the mesh size. We show agreement of the results 

obtained by this microscopic statistical mechanical theory with experimental 

data and the existing scaling arguments. 

PACS numbers: 61.41. +e, 64.60.Cn, 64.60.Kw, 64.70.Pf 
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Polymer gels exhibit a wide spectrum of interesting properties resulting from their me- 
chanical stability, their elastic properties and their ability to swell upon absorbtion of solvent, 
sometimes giving rise to a multitude of swollen phases |1|]. The theoretical approach to such 
systems with cross-link constraints has also been related to the study of the determination 
of protein structures given a set of constraints produced by protein NMR methods 0. 

The statistical mechanics of systems with random cross-links have attracted significant 
attention and have been the subject of extensive debate in the interpetation of the results 
obtained by the different methods f^-Q. In order to take into account the randomness in 
the distribution of the cross-links the replica method was used and the interpretation of 
the multi-replica order parameter was lacking. Recently, Panyukov and Rabin |^ derived 
a physical effective Hamiltonian that allows the calculation of observable parameters such 
as density correlation functions due to randomness in the quenched cross-links and due to 
thermal fluctuations. Besides this seminal contribution, their theory is also consistent with 
the behavior of thermal and quenched fluctuations due to deformation (butterfly patterns) 

EE- 

In the present work we apply the Panyukov- Rabin method to derive the effective Hamil- 
tonian of a system of a blend of two very long homopolymeric chains of different kinds that 



are cross-linked llOHl^l- We consider the case where there is an effective repulsion between 
the two kinds of monomers that leads the uncross- linked blend to segragation at low tem- 
peratures. We assume that the cross-linking takes place at high temperatures where the 
blend is in the mixed state. After the system is cross-linked reduction of temperature will 
not lead to macroscopic segregation because the cross-links hold the pieces of different kinds 
together. Instead a microphase separation will appear. 

The study of this transition so far has been based on a phenomenological Hamiltonian 
introduced by de Gennes [|T0| where the effect of cross-links has been taken into account 
as a Hook-law type of elastic forces that keep the pieces of different kinds together. A 
complete derivation of the effective Hamiltonian starting from a microscopic model, to the 
best of our knowledge was missing. The result obtained here provides a deep understanding 



of the physics of the system and its quantitative predictions are in fine agreement with 
experimental data as discussed below |T^. 

In a microscopic description we consider a Gaussian chain with Kuhn segment a described 
in the Edwards formulation by 



i *U) +— i *i wo -0)1. (1) 

where w^^-* the excluded volume. In our problem we have two kinds of chains, e.g., one white 
chain and one black, that have a tendency for segregation. We denote the position of the 
monomers of one kind by x(i) and the other by y{i). The interaction between the two kinds 
is given by 

n, = -xjdij dj5[x(0-y(j)]. (2) 

The cross-link constraints are introduced by a set of delta-functions and the partition 
function can be written as 

z= Px(z)exp[-7^o(x«)-7^o(y«)-^i(x«,yW)]n^[[^fe)-yOp)]' (3) 

p=l 

where p is a label for a particular cross-link and the number of cross-links in a particular 
realization of disorder. 

Then we use the replica trick to perform the average over all realisations of cross-links. 
However, taking this disorder average it is important to impose an additional constraint: 
that the cross-links are compatible with an existing real conformation @J5J3l- ^o ensure this 
we fix the cross-links to be compatible with a conformation labeled as an additional "zeroth" 
replica. Also with this constraint we ensure that the number of ways to chose N^. cross-links 
should not exceed the total number of available conformations and the conformations over 
which the partition function average is taken are physical. Then the n-th power of the 
partition function averaged over disorder reads: 

(Z") =Jdsfl{J V^-iz) J I)y"(Oexp[-Ho(x"(0) -Ho(y"(0) -Hi(x"(^),y"(0); 

xn^[x"(^p)-y"Op)] > (4) 

p=l ^ 



where the average over all cross-link configurations is given by 

The derivation of the effective Hamiltonian from this model can be done by using the 
grand canonical ensemble with a fixed number of N monomers for each chain. The first 
order parameter that appears in this process is 

rN " 

pi(S)=pi(x°,x\...,x",...,x")= / n'^(^"«-x°), (6) 

where x is a multi-replica 3(n+l)-dimensional probe vector, x"(z) is the position of monomer 
i of replica a, and x° is a probe position in the 3-dimensional space of replica a. The 
corresponding parameter P2(x) is defined for the monomers of the other kind. This parameter 
can be interpreted as the average correlator of conformations of the n replicas with the 
"zeroth" replica which corresponds to an existing conformation at the moment when cross- 
linking takes place. This conformation is also taken to be maximally compact with constant 
real density. 

In order to derive the effective Hamiltonian for our system we need to calculate the 
entropy of conversion from the microscopic to the macroscopic parameter with introduction 
of the fields (pl (x) /2 = pi (x) and 02 (^) /2 = P2 (x) . The part of the effective Hamiltonian 
that does not contain the repuslive interaction between monomers of different kinds Tii, then 
reads 




The chemical potentials pi,p2 are introduced in the grand canonical ensemble formalism to 
fix the number of monomers of each chain to N. They are equal with each other and are 
given by 



where N is the average number of monomers from the one cross-hnk to the next and p the 
mean field density. 

The number of cross-hnks is proportional to which is given by 

Zc = i. (9) 

The excluded volume effect is important for the "zeroth" replica since it fixes the maximal 
number of cross-links at the cross-link saturation threshold w^^^ pN = 1. 

This Hamiltonian and the corresponding multi-replica order parameter do not give ap- 
preciable physical insight and some physical order parameters should be sought. Following 
the methodology of Panyukov and Rabin we perform a procedure for the evaluation of the 
entropy for the conversion into an order parameter related to differnce of the local densities 
of the two kinds. 

Originaly we observe that the mean field values for the densities of the two kinds of 
monomers should be equal because of symmetry. It can be shown that this mean field 
solution is given by a one variable mean field equation 



The single variable s of this differential equation corresponds to the average deviation of 
the n replica conformations from the "zeroth" replica, i.e., the initial conformation at the 
moment of cross-linking. The solution of the mean field equation can be done numerically. 

However, the local fluctuations should account for the differences that give rise to a 
microphase separation pattern. The microphase separation parameter should be calculated 
as a difference between the single replica fluctuation differences. The latter can be calculated 
from 




(10) 



where 




(11) 




(12) 
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If we set 01 (x) = 01 (x) +■01 (x) and 02 (x) = 02 (x) +'02(x), where ■0i(x), ■02(x) the fluctua- 
tions in the multi-rephca fleld, the lowest order estimate with respect to fluctuations for the 
difference between densities of the two kinds is given by 

2 /(n d^^){M^) - M^))'Pmf{s). (13) 

We need to evaluate the entropy for the conversion to the physical single-replica order 
parameters that correspond to the microphase separation in real space. This can be done 
with the introduction of an auxiliary fleld 9 and the n-th power of the partition function 
which reads 

{Z^') = [ Pm«(x«)P^„(x")exp[-xE / c?x-m^(x") + ^ / cix-m«(x")^,(x")" 

x(cxp{-22 / dx0^/(s)[V'i(x) - V2(i)] E^"(x")}) , (14) 

where mo,(x") is the macroscopic microphase separation order parameter as a result of the 
difference between fluctuations of the densities of the two kinds of monomers and ISH is the 
part of the Hamiltonian that depends on these density fluctuations, given by 

A7i[V^i(x),^2(x)] = \j dx[l(V^,2(x) + V^2'(x))]-y[V2V'i(i) + VV2(i)] 

-|[V^?(x)+V'2'(x)]0^^(.)-2ZeV'l(i)V2(i)0^/(«). (15) 

The evaluation of this conversion entropy results into a more physical expression for the 
effective Hamiltonian which reads 

where ruqa is the Fourier transform of the microphase separation order parameter. At this 
stage we see that there is no coupling between replicas so that the replica index becomes 
irrelevant and can be droped. The newly introduced fleld n corresponds to the quenched 
local differences of the densities of the two kinds of monomers associated with the initial 
conformation. The deviation from this random quenched microphase structure raises the 
energy due to the network elastic forces. The form of the correlation function will depend 
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on the thermal correlation function gfq and will be discussed below. The field n defines 
the profile for a reference state with fixed density. This state is very different from an 
ideal reference state without excluded volume, which is collapsed to the dimensions of the 
mesh size and cannot serve as reference to the excluded volume problem problem of an 
incompressible melt. Although an elegant way has been recently proposed for the ideal 
chain with cross-links |]14| this cannot serve as a basis for the solution of the incompressible, 



real melt. 

Thermal fluctuations correspond to probing a particular part of space in a gel with a 
given set of cross-links and corresponding to a time average. The average over different 
disorder realizations, i.e., different sets of cross-links, should be done with the aid of the 
probability distribution of the fleld n which can be calculated according to the principles 
laid by Panyukov and Rabin and will be presented elsewhere since, as it will be discussed 
below, it is not important for the study of the microphase separation transition principles. 

The most important quantity is which can be calculated by a long procedure that 
involves numerical solutions of differential equations and integrals and will be presented 
in detail elsewhere. The result obtained for the large scale limit q ^ is (?(q — 0) = 
O.eOa^iV (p)V2q2. This expession is different from the corresponding correlation function 
for the density fluctuations found by Panyukov and Rabin in that the leading term is pro- 
portional to and not of order unity. 

For short scale behavior, i.e., at scales shorter than the mesh size where the system does 
not feel the cross-links, the density fluctuation correlation function is given by G'o(q) = 
12/(a^q^) for each of the two chains of our system |T^. By setting pi = (1/2) (p + m) and 



P2 = (1/2) (p — m) we can flnd easily the contribution of this term to the free energy, so that 
the effective Hamiltonian to second order becomes 



1.66 2 1/ \ 

-2. . + TTT^ + nix - Xo) 



m. 



2 

q' 



(17) 



.a2Ar (p)V2q2 24 

where (x — Xo) is the shift from the temperature of segregation for the uncross-linked blend. 
The prefactor 1.66 differs from the estimate 36 of de Gennes. The fourth order vertex is 



much more difficult to calculate and is an open question. 

The inverse correlation function of Eq. (17) gives rise to a Hamiltonian of the Brazovskii 
||T6| type A{q ~ qo)"^ + t with a minimum at 

a% = _i/2^ , , (18) 

and the density can be set equal to unity for simplicity. This coincides with the result 
obtained by de Gennes appart from a numerical prefactor. The correlation function diverges 
for q = go aX the transition point r = 0, whereas the quenched fluctuations do not depend 
on temperature, i.e., the proximity to the transition point, so they will not be important for 
the observations near the transition point. 

The prefactor 2.5 should be considered more accurate, because it comes from a micro- 
scopic theory and is in agreement with the experimental value of 2.3 which is found from 
the microphase separation wave length mesurements in a blend of deuterated polystyrene 



with poly(vinyl methyl ether) crosslinked with radiation ITT 



With a statistical mechanical treatment for gels of the sort described here and the en- 
couraging results obtained in this work the way is now open to tackle more chalenging 
problems in the filed of thermodynamics of heteropolymeric gels such as the case which in- 



cludes electrostatic forces and which shows complex behavior and multiple phases |]18[ that 
have remained an open problem in the field for along time. 
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